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ABSTRACT 



Presently, many models of the coronal magnetic field rely on photosplieric 
vector magnetograms but these data have been shown to be problematic as the 
sole boundary information for nonlinear force-free field (NLFFF) extrapolations. 
Magnetic fields in the corona manifest themselves in high-energy images (X-rays 
and EUV) in the shapes of coronal loops, providing an additional constraint that 
at present is not used due to the mathematical complications of incorporating 
such input into numerical models. Projection effects and the limited number of 
usable loops further complicate the use of coronal information. We develop and 
test an algorithm to use images showing coronal loops in the modeling of the solar 
coronal magnetic field. We first fit projected field lines with field lines of constant- 
a force-free fields to approximate the three-dimensional distribution of currents 
in the corona along a sparse set of trajectories. We then apply a Grad-Rubin- 
like iterative technique to obtain a volume-filli ng nonlinear force- free model of the 
magnetic field, modifying method presented in lWheatlandl (120071 ). We thoroughly 
test the technique on known analytical and solar-like model magnetic fields previ- 



ously used for comparing different extrapolation techniques (ISchrijver et al. 



2006 



20081 ) and compare the results with those obtained by presently available methods 



that rely only on the photospheric data. We conclude that we have developed a 
functioning method of modeling the coronal magnetic field by combining the line- 
of-sight component of photospheric magnetic field with information from coronal 
images. Vector magnetograms over the full or partial photospheric boundary of 
the numerical domain could optionally be used. 



1. Introduction 



The ability to build adequate models of the coronal magnetic field is extremely important 
for understanding the physics of the solar corona. The corona is believed to be generally in 
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a force- free (or at least low-/3) state ( lGaryll200l[ ). Destabilization of this state may lead to 
eruptions, with contributing factors including topologi cal properties of the field, such as the 
existence of null points and excessive magnetic twist ( ICanfield et al.lll999l ). The amount of 
energy released in eruptions cannot exceed the amount of free magnetic energy at the time of 
destabilization. Moreover, as the coronal field generally evolves in such a way that its total 
helicity only changes due to helicity flux across the photosphere and into the heliosphere 
( iBerger fc Field! Il984j ). the assessment of helicity at one point in time, such as prior to a 
CME, might be beneficial for studies of the evolution of the corona and heliosphere. Modeling 
of coronal heating is frequently performed as 1-D hydrodynamic (or static) models following 
mag netic field lines using y alues of magnetic field along these field lines as important input 



(e.g., iLundquist et al.ll2008[ ). So this modeling would also benefit from better models for the 
magnetic field. 

The general problem of constructing a force-free magnetic fie ld (hereafter FFF) to model 
the coronal field is formulated as follows (INakagawa et al.lll97ll ). The objective is to find a 
magnetic field B which satisfies the divergence-free condition 

V-B = (1) 

and the force-free equation 

V X B = aB, (2) 

where a is a proportionality constant between the magnetic field and magnetic currenl0. 
Equations ([T]) and ([2]) must be solved for B and a in a volume domain V subject to boundary 
conditions B|gv (or B • n|gv and a|av)- The problem is not in general linear and the solution 
is hence called a "nonlinear force-free field", hereafter NLFFF. Particular cases include a 
linear force- free field (hereafter LFFF) that solves the system assuming a(r) = const, or a 
potential field (where a(r) = 0) which we refer to in the text as Bp. 

Many difficulties arise when solving the problem of constructing a non-linear force-free 
field. The underlying reasons for these difficulties are physical, mathematical and compu- 
tational. Physically, the full vector magnetic field at the lower boundary z = is presently 
obtained only in the photosphere, where plasma forces are significa nt. That is to say. Equa- 
tion ([2D is not appropriate at the lower boundary level (jGaryll200ll ). Also, the component of 
B transverse to the line of sight at the photosphere is subject to an intrinsic 180° ambiguity, 
and measurements of boundary data at the top and side bounda ries of the computational 
domain are not available at this time 
sion of these issues 



see 



Demoulin et al 



19971 . for an extensive discus- 
Typically, assumptions are made about the side boundaries, e.g., a 



^The parameter g has topological meaning associated with the amount of twist in the field, e.g., 
Gold fc Hov3 (Il960l) . 
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field matching a potential source surface model (ISchrijver &: De Rosall2003h is assum ed, and 
there are various methods to resolve the azimuthal ambiguity (iMetcalf et al.l 120061 ) . Math- 
ematically, the system is nonlinear and at the present stage the uniqueness, and even the 
existence, of a solution in general for a given boundary conditions are not proven. Finally, 
there are computational difficulties that have to do with the high instrumental uncertainty 
in the measurements of the transverse horizontal component of the photospheric magnetic 
field and the small spatial scale of current changes, possibly below the instrumental reso- 
lution, in the lower boundary. This uncertainty has more impact than it might seem at 
first sight because B ■ Va = 0, implying a = const along magnetic field line^. Hence, field 
lines must connect points with the same a on positive and negative polarities at the lower 
boundary so the boundaries must have equal amounts of incoming and outgoing magnetic 
flux for each value of a. Noise in a at the lower boundary and limits to the fleld of view pre- 
vent this condition from being satisfied, and the problem is in general ill-posed. Techniques 
e xist for "pre-processing" of t he boundary data to attempt to mitigate this problem (e.g.. 



Wiegelmann et al.ll2006l . 120081 ). 



The existing methods to address the difficulties outlined above do not appear to be 
developed to a level such that photospheric vector magnetograms may be used to reliably 
model the coronal field. Different methods for solving the NLFFF problem, and even different 
implementations of the same method, applied to the same photospheric data, and even the 
same method applied to different polarities of the sa me data, frequently yield results inconsis- 



tent with each other and with the corona l features (ISchrijver et al.ll2006l : iMetcalf et al.ll2008 



Schrijver et al. 2008; DeRosa et al. 20091). Suc h methods are, for ex ample, the magnet ofric- 



van Ba 



tional relaxation (e.g. 
the Grad- Rubin method (e.g.. 



Iegooiierul2004j). optimization (e.g., IWiegelmanru |2004J ) and 



Wheatland 



20071 ). 



Extensive studies are needed to address all of these issues. Hence, a substantial time 
might pass before reliable vector magnetograms consistent with the upper chromosphere 
become available for models of the coronal field. Until they are available another source of 
information is needed for modeling the coronal fields. We propose this source to be coronal 
loops. 

Coronal loops, observed in X-Ray and EUV images, are believed to follow lines of 
the magnetic field, and therefore they should be of help for magnetic extrapolations. Unlike 
vector magnetograms, this information originates in the force-free corona, where Equation ([2]) 
is appropriate. Field lines spread apart with height and so do bundles of coronal loops 



^This follows from Equation ([2]) by taking the divergence of both sides of the equation. 
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(though the field generally expands with height, ind ividual loops are found to have nearly 
constant diameter with height, see iKlimchukI |2000[ ). Consequently, the structure of the 
magnetic field in the corona should be less fine than at the photospheric level so it might 
in principle be better resolved by currently available instruments. Observed loops also give 
an idea about the overall connectivity of the coronal field, which might otherwise be easily 
distorted by even minor noise present in photospheric vector magnetograms and therefore in 
a, as discussed above. 

Even if techniques of processing vector magnetograms are developed to the point that 
NLFFF models are generally reliable, coronal loops as an additional constraint might be of 
great benefit, for example, for studies of energy release in solar flar es. Vector rnagnet ograms 



undergo relatively minor changes during even major flares (e.g., IWang et al.l 120121 . found 
only a fractional change in the transverse component of the field in a small patch of the 
active region during a large X-class flare). On the contrary, the changes in the connectivity 
of the coronal magnetic field can be large-scale and dramatic even in smaller flares. As the 
connectivity of the magnetic field manifests itself in the shapes of coronal loops, the latter 
provide a powerful guide for tracking sudden changes in the field. 

Making use of coronal loops is, however, a non-trivial task. The plasma is optically thin, 
and what is observed by instruments is the integrated emission of all the plasma along the 
line of sight. Extracting individual loops from bundles of overlapping loops is a non-trivial 
image processing task, with the possible exception of isolated loops far aw ay from the core of 



the r egions. Some progress, though, has been made in this direction (e.g.. lAschwanden et al. 



20081 ). Another diffi culty is that all currently existing instruments, with the exception of 
STEREO satellites ( jKaiserl |2005| ). only observe the Sun in one projection, so the three- 
dimensional structure of the loops is not immediately obvious. 

Rece ntly, substantial p rogress has been made in studies of coronal loops as magnetic 
features. Lim et al. ( 20071) fir s t fitte d observed projections of coronal loops with lines of 
a LFFF. iMalanushenko et al.l ( 120091 ) developed a semi-automatic algorithm for such fits, 
applicable to portions of loops, and showed that a values obtained this way statistically 
correlate with a values for a NLFFF model. Progress also has been made in obtaining 



information from the original images. Numerous studies, (e.g., lAschwanden et al.l |2009[ ) 



have demonstrated good results on triangulating loops trajectories using STEREO data. 

Two approaches to the use of image data in magnetic modeling are as follows. Recon- 
structed 3-D loop traj ectories and a values along them may be determined approximately 
using the scheme from IMalanushenko et al.l ( l2009l ). hereafter the MLM09 fit. This provides 
information at least about the 3-D trajectories of some field lines and a in the corona along 
these field lines. Stereoscopically-derived data offers another possibility: the inferred 3-D 
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loop trajectories could be used in conjunction with values of the vector magnetic field at 
the loop foot points. Vector magnetograms are of course prone to the problems outlined 
above. However, in the case of using loops, the field values only need to be accessed at a 
sparse set of locations in the lower boundary. If it is possible to estimate the chromospheric 
magnetic field (assuming, for example, that the field does not change much with height in 
the chromosphere) in at least a few patches in an active region, and provided that these 
patches contain foot points of the stereoscopically determined loops trajectories, then this 
information could be used as in the first approach, but with more accurate results. 

In this paper we propose a new method of constructing a NLFFF using such information 
derived from coronal loops. We also draw attention to the value of coronal loop observations 
for magnetic modeling in general. Such methods might in principle be of use in areas of 
plasma physics other than coronal studies. It might for example be desirable in laboratory 
plasma studies to estimate what kind of a force-free field would have a required topology 
and magnitude of currents. 

The paper is organized as follows. In Section [2] we describe the quasi Grad-Rubin 
scheme enabling us to make use of coronal loops with and without vector magnetograms. In 
Section [3] we discuss various inputs. Section H] describes the general scheme of a set of tests 
of the method and figures of merit obtained. The results of the tests are presented in detail 
in Section [51 Section [6] discusses the results, evaluating how successful the scheme is and its 
value for modeling of the coronal field. 
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2. Description Of The Quasi Grad-Rubin Method 

Suppose there is a domain V with boundary dV and the following are given: 

(a) B • n|av (where fi is the normal to dV); 

(b) a set of trajectories {Ci}fLi in V along which the force-free parameter values {aj}^ 



1=1 



are known (and are constant along each individual trajectory). 



The objective is to find the field B that solves Equation and matches the boundary 
conditions (a) and the volume constraints (b). 



The procedure is iterative and is similar to a Grad-Rubin iteration ( iGrad fc Rubin 



19581 1 ■ It starts with potential field B(°) = Bp as an initial guess for the field and an initial 
guess a^^^ for the force- free parameter, which at each point in the domain is set equal to 
ttj for the closest point in the volume for which an at value is known. Then on every n-th 
iteration the updated cubes B^") and a*^"^ are obtained as follows: 

1. Impose the volume constraints by setting along the trajectories Ci. 

2. Calculate updated field values B*^") from 

V X B(") = a("-^)B("-^) (3) 

subject to the prescribed boundary conditions. This equation is solved using a vector 
potential A*^") such that B*^") = V x A'-"'', so the divergence-free condition is satisfied 
to truncation error. 

3. Calculate an updated set of values for the force-free parameter a^"-*: for every point in 
V, assign a^"'^ = {a^'^~^^) averaged along the field line in B*^") that passes through that 
point. If a field line leaves the domain through any boundary but the lower one, the 
value of a is set to zero along it (in common with the Wheatland 2007 GR scheme). 
This ensures that no currents go off to infinity so that the fields' energy remains finite. 

4. Repeat 1.-3. until B*^"^ ^ b*^""^^ and a^""^ ~ a*-""^) to within a tolerance and therefore 
V X B*^") ^ a'^^^B*^"). 

This sequence is similar to an existing Grad-Rubi n method of solution of the NLFFF 



problem ( IWheatlandl 120071 : IWheatland &: Regnierl |2009| ) . The only difference between these 



two schemes is how the updated a*^"-* cube is calculated in the Step 3. In both schemes, at 
each point in the domain a field line is traced in B'^"^ In the original Grad-Rubin code. 
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positive or negative polarity is picked at the lower boundary; and a^^^ at each point in the 
volume is set to the value at the boundary point in the chosen polarity where it is crossed 
by the field line. The only exception is the boundary itself for all points in the chosen 
polarity where a is kept constant. In the Quasi Grad-Rubin, a^") at each point is assigned 
the average of a from the previous iteration along this field line. The only exception is the 
volume constraint paths where a keeps the volume constraint value. 

Let us consider some particular simple cases. First of all, it is clear that if the initial 
guess for B and a already satisfies V x B = aB, then the scheme keeps it unchanged. If 
the field or a differ from a solution to Equation ([2]) even at one point, the field changes - 
though if it is only one point that is different, convergence is achieved in a single iteration. 
It is also clear that if a = everywhere, currents do not appear; since the scheme averages a 
at every iteration it is incapable of introducing |a| > max(|Q;j|). These cases make sense: if 
the answer is close to the correct answer, convergence is achieved rapidly, and if no currents 
are specified to start with, the scheme does not change the input potential field. 

But what happens in an intermediate situation: currents are known on some, but not all 
fiux tubes in the domain? Can a solution be reached at all? If several solutions are plausible 
given the constraints, which solution (of any) is achieved? There are proofs of existence and 
uniqueness of solution for the force-fr ee problem, a chievable by Grad-Rubin iteration, if a is 



sufficiently small in some sense (see iBineaul Il972l ). However, the range of a is not clearly 



defined and it is unclear whether solar-like fields are within this range. If they are outside 
of this range, that does not mean, the original Grad-Rubin iteration necessarily fails. It is 
unclear if similar proofs exist for the proposed Quasi Grad-Rubin scheme. In the absence 
of this, we consider a process of numerical experimentation to test the scheme. We attempt 
also to determine how many field line trajectories Ci are sufficient to enable convergence. In 
Section O we review some of our experiments. 



3. Different Types of the Input Data for Quasi Grad-Rubin Scheme 

The Quasi Grad-Rubin numerical scheme (hereafter "QGR", in contrast to "GR" for 
Grad-Rubin algorithm) could in principle be used with a constrained at any set of locations 
including the lower boundary. Hence it may be used with vector magnetograms, setting a, 
at the lower boundary to the vector magnetogram derived value. 



a\ 



1 fdBy OB, 



(4) 

2 = 



Bz \ dx dy 

Vector magnetograms may also be used with or without the loop trajectories. We identify 
three different kinds of inputs for QGR: a along loop trajectories; a along loop trajectories 
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and at the lower boundary; and ot at the lower boundary onljo. Traditional schemes are 
designed to work for the last case only. 

If the a values, wherever set, are approximate, this will introduce uncertainties. To 
properly test QGR we try to recover several known fields and use both approximated tra- 
jectories and those drawn from the known field. We note that there may be applications 
of QGR even if it would not work with approximate data, for example for problems of the 
following kind: constructing a magnetic field with twist prescribed along certain trajectories. 
It could also be used with stereoscopically triangulated loops and currents derived from the 
photospheric vector magnetograms, or a values and trajectories obtained by other means. 

In this paper we test various possible inputs and compare the results with the reference 
fields. Table [T] outlines the degrees of freedom available for such tests. We refer to various 
inputs and schemes using this chart, e.g., II. b is QGR applied to volume constraints alone, 
drawn from the reference field. 



Table 1: Possible combinations of different inputs to QGR, in addition to B • n|gv Note that La 
has the same input as the original GR scheme but the algorithm is different, so QGR with La input 
is not equivalent to GR. 





Values of a at 2; = 


Values of a along loop trajectories 


Known 


Unknown 


None 


La 




From the field lines of the model field ("ideal" input) 


Lb 


Il.b 


From the MLM09 approximation derived from 2-D projections 


Lc 


ILc 


of these field lines (realistic coronal input) 



The QGR scheme was implemen ted by modifying an existing GR code (which we refer 
to as CFit version 1.3), described in IWheatland fc Regnierl ( 120091 ). It is a "self-consistent" 
scheme, in that it picks a polarity, finds a force-free solution using boundary data from that 
polarity, hence obtaining values of a everywhere in the volume including the other polarity 
at the boundary; then the a map at the other polarity is updated with the weighted average 
of the values obtained from this new solution and the values that existed before this solution 
was found. It then repeats the cycle using the updated a map from the other polarity. The 
cycles are continued until the two solutions obtained using a values at the opposite polarities 
are consistent with each other to a tolerance. Since QGR does not use the lower boundary 



■^Note that a does not have to be constrained at all points on the lower boundary. 
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in the same way, we discard the switching between cycles and modify the way the a values 
in the volume are calculated (see previous section). Also, instead of using at z = and 
a 2-D array of a values at z = the modified code uses at z = and two 3-D arrays: 
along the trajectories and an initial guess everywhere else and a "mask", i.e. another 3-D 
array with entries either unity along £j or zero everywhere else. These are the two principal 
changes to the CFit code. The calculation of the field on a given iteration, i.e. the solution 
of V X B*^"+^^ = a^'^^'B^"'^ is unchanged. The code uses a vector potential for this step and 
hence the field satisfies V ■ B*^"+^) = 0. 



4. Description Of Metrics for QGR Solutions 



In this paper we try to recover seve r al kn own force-free fields, namely those from 
Schrijver et al.l ( l2006l ) and ISchrijver et al.l ( l2008l ). The iterations are initialized with the 
same potential fields used in the original studies. We also try to construct a NLFFF based 
on a dipole magnetogram and two specified loop trajectories (with a reference field which is 
not known). 

We estimate the quality of the reconstruction and the relative force-freeness of the known 
solutions using metrics of which most have become standard in NLFFF modeling. These are 
as follows. 



• E/E-p versus Ej-cf/Ep: energy in the reconstructed field versus energy of the reference 
field (for a perfect reconstruction the two would be equal). 

• if(B|Bp) versus //(Bjef |Bp): relative hehcitjfl in the reconstructed field versus that 
of the reference field (for a perfect reconstruction the two would be equal). 

• CWsin = (X] I sin^^||J|)/(^ |J|) versus CWsiUref, where |sin6'| = |J x B|/|J||B|: the 
total current-weighted sine of the angle between B and J (for a perfectly force-free 
field this is zero). 

• Metrics of similarity between B and B^ef, normalized to equal unity if B = Bref: 

— Ccs = ;^ ^ [B ■ Brcf/(|B| |Brcf I)] (where N is the number of points in the domain): 
the average cosine of the angle between B and Bref; 

— Cvec = (Yl B ■ Bi.ef)/(^ |B||Bi.ef|): Same as previous but with increased weight in 
regions of stronger field; 



"'By _ff(Bi|B2) we mean the helicity of the field Bi relative to the field B2. 
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E'^ = I - where Em 
between B and Bref; 



^E|B-B 



N 



rcf 



IB 



ref 



the average relative difference 



E'^ = 1 — En, where En = E |B — Brefl/E iBref]: same as previous but with 
increased weight in regions of stronger field. 



We omit metrics for how well V ■ B is satisfied, because the method fin common with 



Wheatland &: Regnierll2009[ ). uses a vector potential to calcul ate the field and hence achieves 
a divergence-free state to truncation error (jPress et al.lll992[ ). 



5. Sample Applications of QGR 

5.1. QGR Solution for a Dipole Field 

The first test case is a simple dipole field aligned in the E-W direction with the North 
half of both magnetic poles having negative twist and the South half of both magnetic poles 
having matching positive twist. This model could be viewed as a simple representation of 
an emerged untwisted flux rope whose foot points became distorted in such a way that the 
field at one (leading) polarity has been inclined more than the field at the second (following) 
one, perhaps d ue to subsurfa ce flows. Such a difference in inclinations is observed for solar 



active regions (lHowardlll99l[ ). 



To construct the field we calculate two constant-a fields confined to half spaces (jChiu fc Hilton 



19771 ) with equal and opposite twist (ao = ±1.57r/L, where L is the size of the domain). 
We draw one field line for £j from each of these and use these field lines (which imitate 
coronal loops) as volume constraints. The initial guess for a is ±ao in the S and N halves, 
respectively. The fields, field lines, and the locations of the constraints are shown in Figured] 

QGR is found to converge to a nonlinear force- free field, as shown in Figure O Achieve- 
ment of the force-free state is shown by the distribution of sin(B, J) peaking at zero. A 
second smaller peak at one is a contribution from the current-free regions of the field, and 
is due to J in these regions being exclusively numeric noise. To support this, the shaded 
histogram on Figure [2] shows the distribution of the same quantity evaluated for potential 
field. The nonlinearity of the field, that is, the presence of different a values on different field 
lines, is illustrated by a distribution of a on a horizontal slice. These additional diagrams 
are shown in Figure |2] (a)-(c). 

The field lines of the solution, interestingly, are found to follow different trajectories 
than the two constraining field lines. In fact, the solution, having a close to the two imposed 



- 11 - 



a = — I.Stt/L o;=1.57t/L Volume constraints 




Fig. 1. — Input data for the dipole field test case: two field lines (thick red) drawn from two 
constant-a fields (left and middle panels). The right panel shows points where the volume 
constraints are applied (i.e., the points through which the trajectories Ci pass). Points in 
the north and south halves of the domain are assigned values of iao, respectively. 

constraints in the two halves, appears much closer to a potential field than the highly twisted 
constant-a fields we drew the loops from! This is shown in Figure |2] (d) which illustrates 
the constraining loops and two field lines of the solution initiated at the midpoints of these 
loops. 
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Fig. 2. — The solution for the dipole test case. The QGR iteration converges as shown in 
panel (a), which displays the change in the free energy of the field consecutive iterations. 
The constructed field B is force-free, as shown in panel (b) by a histogram of | sin(B, J)| for 
Bp (shaded curve) and B (sohd curve). The peak at unity is due to the current-free regions 
in the field, as explained in the previous section. The initial field Bp is current-free and 
B retains current-free regions (in particular where field lines leave the domain through side 
or top boundary). The field B is nonlinear, as seen in panel (c), which shows a horizontal 
slice of a at a height of six pixels in the box (out of 64). In this panel black corresponds 
to negative a and white to positive a. The field lines of B, however, do not match the 
constraints! Field lines of B are shown in panel (d) as solid black while the constraints are 
dashed red lines. Field lines are initiated at midpoints, rather than foot points, of the loops. 
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The fact that the NLFFF solution with currents similar in magnitude as in the lChiu fc Hilton 



( 119771 ) constant-a field appears to have much less twisted field lines is due to at least two 



effects. First, as open field lines are requir ed to carry no cu rrents in the NLFFF model 



constructed by a QGR (in common with the IWheatlandl 120071 . implementation of GR), the 



yolume containing cu rrents has to be contained in the computational domain, while the 



Chiu &: HiltonI ( 119771 ) fields have currents in the entire half space. As the volume containing 
currents in the NLFFF is smaller, the current density has to be stronger for the field lines to 
have similar shape. Secondly, currents which run in opposing directions close to each other 
might counteract each other in infiuencing the shape of field lines. 

To illustrate these effects we perform a simple numeric experiment. We repeat the 
computation but set = ±/ao on the constraining paths for several values of /. The 
results are shown in Figure [3], first column. We also repeat the experiment for a box of 
smaller size (cropped on the sides and the top) but with otherwise identical setup (Figure [31 
second column). Finally we repeat the computation with the same lower boundary but 
different volume constraints: a field line in the lower half of the domain is drawn from a 
constant-a field with the same sign of a as the upper half but smaller in magnitude, with 
a = — ao/4 = —3n/8L (Figure El third column). For the original setup, the best match 
between the constraining loops (red dashed curves) and field lines of the solution (black 
solid curves) initiated at the mid points of these loops is achieved with f ^ 10 — 12, while if 
the same (or at least similar) field lines are required to exist in a field with currents confined 
to a smaller domain, the best match between the loops and the solution is achieved with 
/ ^ 14. If both loops have a of the same sign, the best match is achieved for / ^ 6 — 10. (No 
steady solution is found for / = 14 for the setup in the first and third columns; in this case 
the QGR solution continues to oscillate. The solutions for / = 14 in the second column and 
/ = 12 in the third columns exhibits oscillations but of small magnitude. Later we discuss 
these oscillations which may result from the input being inconsistent with a NLFFF solution 
and describe a procedure to damp them. In this section, however, our point is to illustrate 
the significance of the scale factor /.) 

The scale factor / cannot be evaluated a priori, but it may be estimated using observ- 
ables, i.e., coronal loops. This can be done by minimizing the difference between the shapes 
of coronal loops and field lines of different solutions corresponding to different values of /. 
In the next two sections we demonstrate that scaling factors obtained this way are indeed 
a proportionality constant between a on lines from a NLFFF and a of the approximation 



of the se field lines by lines of constant-a fields of the type constructed by I Chiu fc Hilton 



(Il977f ). 
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Fig. 3. — Experiments on the dipolar test case to illustrate the effect of opposing currents 



and domain size on the scaling factor /. See Section [5?T] for description. The notation is the 
same as in Figure [2j 
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5.2. QGR Solutions for Low & Lou Fields 



In this section we try to reconstruct the test field from ISchrijver et al.l (j2006[ ) (the 



"reference field" further in this s ection) using QGR . This test case is a member of family of 



analytic NLFFF's introduced bv ILow fc Loul fllQQOh . 



For the first QGR test we use field lines of the reference field as trajectories Ci and take 
ai values corresponding to the correct a values for the reference field. Physically, this could 
correspond to stereoscopically derived loops with chromospheric vector field known around 
their foot points. These data are hard to obtain at present; we use them mainly to test QGR 
alone, that is, on the "ideal" data not contaminated by measurement errors. We use 113 
randomly selected field lines £i,ref (out of a bigger sample, chosen to be closed field fines, i.e. 
with both foot points on the lower boundary). We calculate ajref numerically everywhere 
in Bref and evaluate it along these 113 field lines. We discard all but 27 of these, retaining 
the ones that may be well fitted using MLM09, to make this test consistent with the next, 
realistic test presented later in this section. We construct two solutions of the same size as 
the reference field (64^ pixels), with and without the additional constraint of vector field 
data at the lower boundary (i.e., scheme s Lb and II. b from Ta ble [H). Figures of merit are 



calculated in the same subdomain as in ISchrijver et al.l (120061 ). This center subdomain is 



also used to estimate a best-matching scaling factor /. 

In the second test we determine £j and ctj based on a fit to the resulting MLM09 field 



derived from the normal field at the lower boundary as in iMalanushenko et al.l (120091 ) . We 
calculate QGR solutions using these data with and without vector field data at the lower 
boundary (I.c and II. c from Tabled]). This tests the applicability of QGR to realistically 
available coronal data. It is not obvious that using approximations to loop trajectories and 
approximate values of a along them is sufficient to create a field model at least as good as 
those derived from vector magnetogram data. To investigate this we project the same 113 
field lines as for the ideal data to the ^ = plane to simulate the appearance of loops in the 
plane of the sky. We treat these 2-D projections as synthetic loops and use them to obtain 
a values aj^MLMog along trajectories A^mlmoq using the MLM09 fit procedure. We discard 
loops which upon visual examination are poorly fitted, which leaves 27 loops that appear to 
have a near perfect fit. The trajectories of these loops and their MLM09 approximations are 
shown in Figure |H 



Keeping in mind the results from Section 15.11 we calculate several solutions for II. c 
with seven different scaling factors / applied to input a but with the calculations otherwise 
identical. Three of these solutions (corresponding to / = 1.23, 1.69 and 2.15) are shown in 
Figure [5l top row. For each of the solutions we estimate how closely the field lines match the 
synthetic 2-D loops used to construct the trajectories Ci. We calculate the average distance 
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Fig. 4. — Results of MLM09 applied to the reference field from lSchrijver et al.l (120061 ) . used 
for I.c and II. c inputs as trajectories Bref field lines (red) and MLM09 field lines (blue). 



between the projected loops and the corresponding field lines of the solution. The result 
(as a function of /) for all seven solutions is shown in Figure [5] (bottom left plot). The 
solution for / ^ 1.69 is the closest match to the loops. The bottom right panel in the same 
figure is a scatter plot of a^t and ttref for individual loops. The fit appears to underestimate 
the value of a. The underestimation factor is remarkably close to 1.69 (dashed line on the 
same plot). This suggests that such underestimation could be derived a posteriori from the 
observed loops. 

The results of our tests (I.c and II. c with / = 1.69) are summarized in Table [2] and 
Figure |6l Figure [7] shows that convergence is achieved for all solutions. We conclude that 
QGR is able to reconstruct the reference field, including the shape of field line^, the structure 
of the currents and the distribution of a, remarkably well. The figures of merit show that 
the QGR reconstructions are at least as close (and typically closer) to the right answer as 
other methods. In particular, the smallest estimate for the free energy we obtain is closer 
to the correct ans wer than any o f the e stimates based on the vector field boundary values 
alone reported by ISchrijver et al.l (120061 ) . 



^Excepting the field lines leaving the computational domain, which are required to carry no current in 
our scheme, but in fact do carry currents in the Low & Lou field. 
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Fig. 5. — Top row: several QGR solutions for the reference field from lSchrijver et al.l ( 120061 ) 
for input II. c (see Table [1]) corresponding to the same input but different scaling factor / for 
a. Red dashed lines show lines of B^f projected onto z = Q and used as loops which were 
approximated by MLM09 to construct trajectories £j. Black lines show corresponding lines 
of B (initiated at midpoints of the lines of Bref). Field lines in the core of the domain are the 
most affected by the choice of /. Bottom left: average distance between loops (projections 
of lines of Bref) and corresponding lines of B in projection onto the plane of the sky {z = Q). 
The difference is smallest for / ^ 1.69. This coefficient is remarkably close to the scaling 
coefficient between Oref and agt from MLM09, shown at the lower right panel as a dashed 
line. Diamonds show a^ei of the individual field lines of Brcf versus agt of the MLM09 
approximations of these field lines. 
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c 



cs 



E„ 



Er, 



CWsin 



E/Ef 



Reference field 
1.00 



1.00 



1.00 



1.00 



0.01 



1.24 



1.00 



Quasi Grad-Ruhin with vector magnetograms 



La 0.99 


0.93 0.80 


0.64 


0.03 


1.27 


0.67 


Lb 0.99 


0.96 0.83 


0.70 


0.02 


1.19 


0.81 


I.c 1.00 


0.97 0.88 


0.75 


0.02 


1.23 


0.91 


Quasi Grad-Ruhin 


with loop trajectories alone 










ILb 0.99 


0.96 0.81 


0.68 


0.02 


1.17 


0.77 


ILc 0.99 


0.97 0.86 


0.78 


0.02 


1.23 


0.93 


Ranges reported in 


Schriiver et al. (2006) 










0.94-1.00 


0.54-0.91 0.48-0.92 


-2.2 - 0.66 


0.03- 


-0.57 0.82- 1.14 




Potential field 












0.86 


0.87 0.50 


0.44 




1.00 


0.00 



Table 2: Metrics (defined in Section Hj) fo r diffe rent QGR solutions with different types of 
input data applied to the ISchrijver et al.l (120061 ) test case. The values for I.c and II.c are 
reported for the optimal solution with / = 1.69. The values for B^ef and Bp are shown for 
comparison, an d so a re the ranges of values for different NLFFF extrapolations reported in 
Schrijver et al.l ( l2006l ). Relative helicity is stated in fractions of that of the reference field. 
For notation, refer to Table [T] 
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Fig. 6. — Reconstruction of Low & Lou field from lScfirijver et al.l ( 120081 ) using schemes II. b 
and II. c, i.e. QGR with ideal and realistic loop input, that is, reconstructed from 2D loop 
projections (see Tabled]). Panels (a)-(i): field lines, hne-of-sight integrated magnitude of 
current and horizontal slices of a for Bref and B for II. b and II. c. Panel (j): field lines of 

"R„ rail ^cAA 1 
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Fig. 7. — Energy E/Ep at each iteration as demonstration of the convergence of the QGR 
iteration for the Low & Lou field test case. The energy on these plots is shown for the entire 
domain (while Table [2] reports the numbers for the middle sub-domain identical to the one in 



Schrijver et al.ll2006l . as do the other tables in the manuscript). Top row, from left to right: 



La, Lb and Lc solutions (refer to Table [H for notation). Bottom row, from left to right: II. b, 
II.c. 
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5.3. QGR Applied to a Solar-Like Field 



In this section, we investigate wlietlier QGR is applicable to solar data. The lLow fc Lou 



(119901 ) family of fields has axial symmetry which is not in general observed in active regions, 
and both magnetic field and current vary unrealistically smoothly through the lower bound- 
ary by comparison with vector magnetogram data. Hence, we repeat the experiments from 
the previous section, but with choos ing a more rea li stic s olar-like field as Bj-ef. We use 



two particular NLFFF solutions from ISchrijver et al.l (120081 ) . who presented NLFFF recon- 



structions of the coronal field for AR 10930 before and after a major flare using several 
extrapolation methods applied to Hinode vector magnetograms. They found that the ex- 
tra polations wh i ch be st matched observed coronal features were GR solutions obtained with 



the IWheatlandl (120071 ) code, using a values from the positive polarity of the magnetograms 
(hereafter Wh+p). We use those as our reference field. These solutions also had the largest 
free energy of all extrapolations. Another advantage of these fields for our study is that they 
use the same boundary conditions and nearly the same numeric implementation as the QGR 
scheme. We emphasize that the objective is in this case not to create a realistic representa- 
tion of coronal field but to test the new algorithm on a known NLFFFs that are expected 
to more closely resemble the coronal field overlying a solar active region. 

For both pre- and post-fiare reference fields we select random sets of field lines and 
evaluate (a) on each field line. These field lines are used as trajectories £j in lib set-up and 
their projections onto z = plane are used as loops for II. c set-up. We als o use at the 



l ower boundary and start with the same initial Bp as the other schemes in ISchrijver et al. 



(120081 ). We perform final tests on the full-sized domain but determine the scaling factor / 
on the domain down-sampled by a factor of 0.5 (this is done to speed up computations and 
to allow the possibility that currents have structure finer than the grid size, which is likely 
to be the case for real data). 

For both sets of loops (volume constraints for the down-sampled domain were only 
imposed at a small fraction of pixels, at about 1.7% of the current-carrying volume and the 
bigger one covered about 6.6% of the current-carrying volume) the solution for the pre-fiare 
configuration does not converge with the II. b inputs. Instead it enters a remarkably stable 
oscillatory cycle with a period of ~ 50 iterations. This cycle develops at ^ 200 iterations, as 
shown in Figure [HI We ran the code for a few thousands iterations to verify that the cycle 
is indeed stable. As energy slowly increases, a sheared arcade forms similar to the one in 
Bref; but as the energy reaches its maximum and becomes most similar to Bref, the field 
experiences drastic changes. Some of the current-carrying field lines rapidly "escape" the 
domain via the y = —100 boundary, what changes the a values on these field lines {a is set 
to zero), as explained in Section [2j In the stage of the cycle with the lowest energy most of 
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the field lines from the core of the region connect to the y = —100 boundary and so carry 
no currents. This may be a valid force-free solution, though it is not consistent with the 
volume constraints which require a to be non-zero in some points in the volume. When the 
a values from the volume constraints are reimposed again at each iteration, the currents 
gradually build up again and the cycle repeats. The escape of field lines does not represent a 
physical evolution of the field as the iterations are not related to any physically meaningful 
time-like variable. The same oscillatory behavior is found in numerous experiments with 



this particular test-case. ISchrijver et al.l (120081 ) report that the pre-flare solution which we 



use as Bref did not fully converge either; it kept oscillating. 

Below we discuss factors possibly causing the oscillations and a way to damp them. 
These factors are: (1) numerical noise in a and therefore in the electric currents that appear 
even in current-free areas and (2) deviations of the input data from a force-free field (as we 
discuss below, Bref in this case is not exactly force- free even at full resolution). The damping 
that we consider allows the calculated field to have small variations in a along field lines as 
well as small deviations of the solution from the volume constraints. 
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Fig. 8. — Values of CW sin and E / Ep in the center of the domain (the same region as used in 



Schrijver et aLl ( 120081 )) for the QGR calculation for one of the solar- like fields demonstrating 



oscillatory behavior. Different stages of this cycle are discussed in the text. The values of 
i^ref and CWsiUref are shown as dashed lines. 
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The first factor is the influence of numerical noise when solving V x B*^"+^'' = a ("^B 
(Step 2 in the algorithm in Section [2]), especially around sharp edges in a^'^\ and with the 
artifacts introduced the by Fourier transforms around these edges. These effects introduce 
noise in the a*^""^^-* values obtained in the next step. Figures 19) and ITO] clarify the amount of 
such noise and its relative size to the signal. In areas of closed field larcf] ~ 0.8 arcsec^^ and 
in the areas with open field (and hence no currents) \aj-ei\ ~ 5 x 10'^ arcsec^^. The fiux- 
weighted distribution of a evaluated numerically in the current-free region has half width at 
half maximum of ~ 10^^ arcsec"^. 

Horizontal slice of a,,, at z = 6.'6 arosec Profiles of a^^, and a ^ at x = 0.0 arcsec, z = 3.8 arcsec 




Fig. 9. — Left panel: image of horizontal slice of a^ef in the pre-flare Wh+p close to the lower 
boundary (the grayscale goes from —0.8 to 0.8 arcsec"^) and two profiles of Oref (solid line 
on both profiles) and Opot (dashed line on both profiles) in this slice. Top right panel: areas 
with significant currents. Bottom right panel: areas with no currents in both Bj-ef and Bp. 
Variations in a due to numerical uncertainties are |a| ^ 0.005 arcsec"^. As we use the same 
numeric solver, the noise in our case is expected to be of the similar nature and magnitude. 



The second factor is errors in the volume constraint data. In the case discussed in 
this section, Bj-ef has significant non-zero magnetic forces from the perspective of QGR. 
As discussed in Section [2l for convergence QGR requires not only that the Lorentz force 
is small everywhere in the volume, but also that that the integrals of the Lorentz force 
along the field lines are small. These conditions are not met for Bref: as shown in Figure fTTl 
a changes substantially along central field lines. This is due to the Wh+p solutions themselves 
not converging preci sely during the GR iteration used to calculate them (as mentioned in 



Schriiver et all 120081 ). 
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Fig. 10. — Histograms of Oref evaluated numerically on closed field (black line) and open 
field regions (gray line) in the pre- flare Whpp. No currents are allowed on the open field 
by the GR scheme, used to calculate Bj-ef. Hence, the gray line shows the numerical noise. 
The distribution of the noise, evaluated from this plot, has a half width at half maximum of 
~ 10~^ arcsec""*^. 

To account for these issues we introduce two uncertainty thresholds: Actgrr and Aonoise- 
The first allows the solution to have a values slightly different from the a values imposed 
along loop trajectories and the second allows small variations of a along field lines to damp 
numerical noise. A revised algorithm is formulated as follows, with the modifications relative 
to Section |2] in bold. 

1. Impose the volume constraints by setting along loop trajectories £j, but 
only at points satisfying 

2. Calculate updated field values B^") from Equat ion |3] subject to the prescribed boundary 
conditions. This equation is solved using a vector potential A*^"'^ such that B*^"^ = 

V X A("), so the divergence-free condition is satisfied to truncation error. 

3. Calculate an updated set of values for the force-free parameter a^"^: for every point 
in V, assign a^^^ = averaged along the field line in B*^") that passes through 
that point, but only at points satisfying > Aanoise- Otherwise 
retain the value of a from the previous iteration. If a field line leaves the 
domain through any boundary but the lower one, the value of a is set to zero along it 
(in common with the Wheatland 2007 GR scheme). This ensures that no currents go 
off to infinity so that the fields' energy remains finite. 

4. Repeat 1.-3. until B*^") ^ b*^"^-'^^ and a*-"-* ~ to within a tolerance and therefore 

V X B(") ^ a(")B("). 
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Fig. 11. — The variation in a in the pre-flare Wh+p field used as a solar-hke test case. Left 
column: selected field lines of Bref (dashed red) and stream lines of Jref initiated at points 
along these field lines. Right column: profiles of ctref along these field lines. These panels indi- 
cate both small-scale and large-scale variation in ctref significantly above the noise threshold. 
The values of {a^ef) and (aref) ± o" are shown as dashed and dotted lines respectively. 

For the case of the Wh+p field we choose Acterr = Aanoise ~ 5 x 10~^ arcsec"^. For 
comparison, significant a for the pre-fiare Bj-ef are |aref I 0.8 shown in Figure [H 

Two exceptions are the II. b pre-fiare solution on the downsampled domain for fewer loops 
case and II. c pre-fiare solution on the full-size domain for more loops case, for which the 
error threshold is increased to 5 x 10~^ arcsec"^ in order to damp oscillations. Figure [17] 
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shows convergence plots for the II. b and II. c cases for the case with the fewer loops. 

The correction factors / are determined in the same way as in the Section 15.21 We find 
a best-fitter / = 4.0 for both pre- and post-flare data. In both cases this factor matches the 
coefficient between aj-ef and agt (see Figures [12] and US]) for individual loops, which provides 
further evidence in favor of this method of estimating /. 
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Fig. 12. — Left panel: average distance between loops from the pre- flare Whpp field and 
corresponding lines of B for the QGR solution in the same manner as in Figure [S] The 
factor / = 4.0, which yields the best-matching solution, is again close to scaling coefficient 
between agt from MLM09 underestimates ttref- Right panel: a scatter plot of Oj-ei and agt 
for individual loops and a line with the slope which equals to the best-matching /. Vertical 
error bars indicate the average variation of a along loops in the reference field. 

The results for the pre- and post-fiare reference fields are summarized in Tables |3lll] and 
Figures [TllfTBl In each case, the QGR reproduces the overall shape of field lines and the 
large-scale features of the current distribution. The reconstructed fields have > 50% of the 
free energy and > 25% of the relative helicity of the reference fields. The reconstructions 
using half resolution data are slightly inferior to those using full resolution data, based on 
the metrics in Tables |3]|11 but the solutions still reproduce at least half of the free energy 
and the quarter of helicity of the reference field and the large-scale structure of currents. 
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Fig. 13. — Same as Figure [121 but for the QGR solution for the post-flare reference field. 
The best-matching scaling factor / is found to be the same as for the pre-flare field. 



Discussion and Conclusions 



In this study we demonstrate that coronal loops provide a useful source of information for 
determining the structure of the coronal magnetic field. While the observed loops do not cover all 
of the coronal volume, they provide information about the shape of the coronal field lines, which 
boundary data alone lack. 

We demonstrate a method that constructs nonlinear force- free fields using line-of-sight magne- 
tograms and coronal loops observed in the plane-of-sky projection. This may mitig ate the problems 



NLF FF schemes encounter with currents determined from vector magnetograms ()Demoulin et al 



19971 ). The loops are first approximated by lines of constant- g fields, with differen t a va lues for 
each loop. This is done using an existing scheme developed by iMalanushenko et al.l (120091 ) which 
we refer to as the MLM09 fit in this paper. The approximate a values along the approximate loop 
trajectories ar e treated as volume constra ints in a quasi Grad-Rubin algorithm, using the code 
modified from IWheatland fc Regnieii (j2009l ) . 



The method, which we refer to as the Quasi Grad-Rubin method (or QGR) is tested on several 
nonlinear force-free fields and the results demonstrate good performance of the method. While 
traditional extrapolations of coronal magnetic fi elds have been foun d to provide poor matches to 
coronal features observed in X-rays and EUV (jPeRosa et al.l l2009l ) . the fields created by QGR 
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are constructed with the effort to match observed coronal features and thus may provide a more 
realistic model of the actual coronal magnetic field. 

The problem of constructing a nonlinear force-free field is typically viewed as a boundary 
value problem requiring an extrapolation of the field from the boundaries to the volume of the 
corona. However, thorough this paper we purposely avoid referring to QGR as an extrapolation 
scheme, because it is not. It is a mixture of extrapolation of magnetic field and interpolation of 
electric currents. The reason we tend to view the step of filling the volume with a values as an 
interpolation-like procedure is as follows. At each iteration, a is averaged along lines of the field 
at the present iteration. If the solution has not converged yet, this field is different from the one 
on the previous iteration, so a is averaged across lines of the field of the previous iteration (see 
Section [2]) . Substantial differences in a values on such field lines result in extreme values getting 
"spread" through the volume and the values of a being "smoothed" along field lines. Smaller 
differences between the fields from two consecutive iterations should result in a smoothing out on 
shorter distances. So on each consecutive iteration n the a values are smoothed across field lines 
to a distance which depends on the magnitude of the angle between B*^"') and B("~^)||J("\ and 
therefore on the Lorentz force at the n-th iteration. The process therefore results in a smooth 
distribution of a in the volume and decreasing the Lorentz forces implies smaller-scale changes in 
this already smooth distribution. The described scheme cannot produce a bigger in magnitude than 
the volume constraints and it tends to produce a smooth transition of a between these constraints 
in such a way as to minimize Lorentz forces. This explains the interpolation-like nature of QGR 
with respect to a. This scheme does not resolve fine structure of the "interpolated" variable (in 
this case, a), as no interpolation scheme can, but it successfully approximates general trends, as 
expected from an interpolation scheme. 

We also develop a way to deal with the uncertainties in the input data and the numeric noise. 
The uncertainties in the observables produce inconsistency with a force-free solution. As such 
uncertainties are expected, this is an important feature of the method. QGR in the form described 
in Section 15.31 allows the volume constraints not to be re-imposed if the average a on a given field 
line which passes through a given constraint point is within a small prescribed amount of of 
the constraint. This means that a magnetic field which is force- free but imperfectly matches given 
volume constraints for a would not be changed by the method. As any numerical scheme, QGR is 
also prone to numerical noise. We are able to determine the range of this noise for a given problem. 
The method assumes that a below this noise level is numerically not distinguished from zero. It 
also assumes that average a along a given field line may vary within this noise range. It therefore 
does not replace a by the newly defined average along the field line if that average differs to less 
than the numerical noise threshold from the previously determined value. 



We noted that fitting loops with lines of IChiu &: HiltonI (jl977l ) constant-a fields results in 



the underestimation of a and verify that it is at least partly due to a difference in the size of the 
volumes which contain currents (finite in reference cases and half space in the fields used for fitting, 
see Section 15. ip . We determine that the underestimation coefficient is roughly the same for most 
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loops and that this coefficient may be determined from observables (projected loops). Applying 
this determined coefficient leads to a good match between the reference field and the model, as 
demonstrated in Section [H A rigorous proof of the nature of such a coefficient and its analytic 
evaluation are subjects of future studies. 

We do not find substantial differences when reconstructing a reference field on a full size domain 
or on a down-sampled domain. This might be due to the smoothness of fields created by the QGR 
due to its interpolation-like nature for the reasons discussed above. This is an important result, as 
the test case in Section 15.31 has structure of currents finer than the grid size in the down-sampled 
test, which may also be the case when modeling coronal fields. 

While developed for currents approximated from loops in EUV and X-ray images, QGR yields 
better results when currents are measured exactly, e.g., from vector magnetograms (II. b inputs in 
Sections 15.21 and 15. 3p . This gives hope that as vector magnetograms become more applicable for 
NLFFF modeling (at least in the cores of active regions), the performance on the II. b level could be 
achieved. This method could also benefit from the exact knowledge of the 3D shape of the loops, 
e.g., drawn from STEREO and SDO satellites combined. 

Overall we find that the method developed in this paper is able to recover over half of the 
free energy and over a quarter of the helici ty for the solar-like test case fields, which is more than 



was r eported for previously tested methods (IMetcalf et al.ll2008l : ISchriiver et al.ll2008l : iDeRosa et al 



20091 ). The method recovers large-scale features of the field well, such as structure of currents, shape 
of field lines and the connectivity of the field, but it fails to resolve fine structure. We nonetheless 
find that the large structure determines at least half of the free energy and a quarter of the relative 
helicity and therefore QGR may be used to provide estimates of these quantities. 
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Cec Ccs I -Em CWsin E / i/(B|Bp) 

Half resolution 

Reference field 

1.00 1.00 1.00 1.00 0.35 1.31 1.00 



QGR with loop trajectories alone, fewer loops case 

lib 0.98 0.98 0.83 0.84 0.37 1.16 0.62 

II.c 0.97 0.97 0.79 0.80 0.32 1.20 0.36 

QGR with loop trajectories alone, more loops case 

Il.b 0.98 0.98 0.82 0.83 0.34 1.23 0.60 

II.c 0.97 0.97 0.77 0.76 0.30 1.27 0.21 

Potential field 

0.86 0.94 0.62 0.70 — 1.00 0.00 



Full resolution 

Reference field 

1.00 1.00 1.00 1.00 0.24 1.32 1.00 

QGR with loop trajectories alone, fewer loops case 

Il.b 0.98 0.99 0.85 0.86 0.11 1.18 0.64 

II.c 0.98 0.98 0.80 0.81 0.07 1.27 0.43 

QGR with loop trajectories alone, more loops case 

lib 0.98 0.99 0.83 0.84 0.09 1.26 0.62 

II.c 0.97 0.97 0.77 0.77 0.08 1.30 0.29 

Potential field 

0.86 0.94 0.62 0.70 — 1.00 0.00 

Table 3: Metrics for the pre-fiare reference field. The numbers for II.c solution are reported for 
the / = 4.0. The downsampled fewer loops Il.b case and full resolution more loops II.c case are 
unstable for AcKerr = 5 x 10~^ arcsec"^; the reported values in these cases are for AcKerr = 5 x 10"^ 
arcsec"^. 
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Fig. 14. — QGR solutions for the pre-flare Whpp reference field in full resolution for fewer 
loops case using schemes II. b and II. c (QGR with ideal and realistic loop input, that is, 
reconstructed from 2D loop projections — refer to Figured]). Panels (a)-(i): field lines, 
line-of-sight integrated magnitude of current and horizontal slice of a for Bref and B for II. b 

on Tnoc! T-l-r^ ( oil ■fiol/^ iTnoo ot-o +t*o/^o/^ ■fyr^rm 
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Field lines of B^gf Ij^gj a^„f at z= 2 pix 
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Fig. 15. — QGR solutions for the preflare Whpp reference field in half resolution for more 
loops using schemes II. b and II. c (QGR with ideal and realistic loop input, that is, recon- 
structed from 2D loop projections — refer to Figured]). Panels (a)-(i): field lines, line-of- 
sight integrated magnitude of current and horizontal slice of a for Bref and B for II. b and 

R„ /'oil fiolr^ 1 
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Cvec Ccs I -Em CWsin E / i/(B|Bp) 



Half resolution 

Reference field 

1.00 1.00 1.00 1.00 0.13 1.16 1.00 

QGR with loop trajectories alone 

Il.b 0.99 0.99 0.86 0.87 0.10 1.07 0.48 

II.c 0.99 0.99 0.88 0.87 0.07 1.13 0.63 

Potential field 

0.94 0.97 0.76 0.80 — 1.00 0.00 



Full resolution 

Reference field 

1.00 1.00 1.00 1.00 0.17 1.14 1.00 

QGR with loop trajectories alone 

Il.b 0.99 0.99 0.89 0.88 0.13 1.09 0.50 

II.c 0.99 0.99 0.89 0.88 0.10 1.14 0.69 

Potential field 

0.93 0.97 0.75 0.80 — 1.00 0.00 



Table 4: Metrics for the QGR results for the post-flare reference fleld. The numbers for II.c 
are reported for the / = 4.0 solution. For notation, refer to Table [H 
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Fig. 16. — QGR solutions for the postflare Wh+p reference field in full resolution using 
schemes II. b and II. c (QGR with ideal and realistic loop input, that is, reconstructed from 
2D loop projections — refer to Figure [1]). Panels (a)-(i): field lines, line-of-sight integrated 
magnitude of current and horizontal slice of a for Bj-ef and B for II. b and II. c. Panel (j): 
field lines of Bp (all field lines are traced from the same starting points). Panels (h), (I): 
line-of-sight integrated volume constraints for II. b and II. c. 
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Fig. 17. — Energy E/Ep at each iteration as demonstration of convergence of different 
schemes for the preflare and postfiare Wh^ field fields. These particular plots correspond 
to the downsampled datacube for the case with fewer loops. 



